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We present a detailed numerical investigation of the tunability of a diffusive random laser when Mie resonances 
are excited. We solve a multimode diffusion model and calculate multiple light scattering in presence of optical 
gain which includes dispersion in both scattering and gain, without any assumptions about the (3 parameter. 
This allows us to investigate a realistic photonic glass made of latex spheres and rhodamine and to quantify 
both the lasing wavelength tunability range and the lasing threshold. Beyond what is expected by diffusive 
monochromatic models, the highest threshold is found when the competition between the lasing modes is 
strongest and not when the lasing wavelength is furthest from the maximum of the gain curve. 


Random lasers (RL) are mirror-less lasing systems 
which have attracted a lot of interest due to their struc¬ 
tural simplicity. Nowadays they have been studied in 
a vast variety of scattering systems ranging from semi¬ 
conductor powder to biological tissue and biocompat¬ 
ible materials [1]. Random lasing originates from a 
complex out-of-equilibrium phenomenon with rich mul¬ 
timodes features [2] and surprising statistical features 

[3] . Despite its potential for practical applications [1], 
random lasing technology is still in its infancy with pi¬ 
oneering applications such as low coherence light source 

[4] and biosensing [5] . One of the factors that has limited 
practical applications is the difficulty of controlling the 
frequency and directionality of the emission. In conven¬ 
tional lasers the lasing emission can be tuned by engi¬ 
neering the high finesse cavity which provides the feed¬ 
back and thus defines the lasing mode. Instead, feedback 
in RL is provided by multiple scattering and the lasing 
emission properties are determined by the complex inter¬ 
play between gain and losses. Recent experiments have 
shown lasing emission controlled by exploiting scattering 
dispersion via resonant scattering sustained by spherical 
particles [6, 7] or by gain dispersion achieved by artifi¬ 
cially increasing absorption in a spectral band [8]. Active 
tuning of the lasing properties has also been achieved by 
shaping the pump profile to selectively excite one or a 
few lasing modes [9-11]. 

Different theoretical approaches to model random las¬ 
ing action have been developed which combine multiple 
scattering and gain. For uncorrelated random systems 
in which interference between the scattered waves can 
be neglected, diffusive models are very accurate even in 
presence of optical gain [12, 13] and they provide the 
time evolution of the lasing process and a smooth lasing 
spectrum with no spiking lasing behaviour [14]. The ra¬ 
diative transport model with gain can also be solved for 
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instance with Monte Carlo simulations which consider 
a random walk of photons [15, 16] and in which am¬ 
plification of single paths can be important in defining 
the spectral properties [17], and by solving the complete 
radiative transfer equations [18]. These approaches al¬ 
low the study of large systems (> 100s mean free paths) 
and geometries similar to real experiments. Recently, 
more complete models including field calculations and 
interference effects have been developed, based for exam¬ 
ple on Maxwell-Bloch equations [19, 20] and the finite- 
difference time-domain solution of Maxwell equations in 
nonlinear media [21-23], but limited to volumes of a few 
wavelengths cubed. 

In this letter, we investigate the range of tunability 
of a diffusive RL with intensity feedback. We solve the 
diffusion equations for typical experimental configura¬ 
tions by including the full spectral dispersion of both 
scattering and gain, beyond the stationary case [12] and 
few modes model [14], and we calculate the full spatial, 
spectral and lasing dynamics, including mode competi¬ 
tion and threshold variation. 

The relevant quantities in a diffusive model are the 
transport mean free path £ t , the diffusion constant D = 
£ t v/ 3, and v the speed of energy in the medium, which 
can all be modulated by resonant scattering [2 ]. The 
optical gain is provided by organic molecules, rhodamine 
6G in this case, described by the stimulated emission 
cross-section cr e , the absorption cross section er a , and 
the lifetime of the excited state r. The molecules com¬ 
posing the gain can be brought to the excited state by 
an intense pump laser. Here we label N\ the density of 
molecules in the excited state and N the total molecule 
density. We model a translational invariant slab geom¬ 
etry similar to most experiments with the ID diffusion 
equation for the pump W p (x, t) and the emitted light 
W' L {x,t). These equations are coupled to a four level 
system gain described by one rate equation for the ra- 
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diative transition. The set of equation is the following: 
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where the emitted light is discretised in spectral bands 
A i, and i = 1 ..n. I p (x,t) is the intensity temporal pro¬ 
file of the pump, 4 is its extinction length and v is the 
velocity of light in inside the medium. The spontaneous 
emission (fluorescence) spectrum defines the quantities 
(j) 1 , with Y2 4> l = 1- When compared to previous works 
[8, 13], the spontaneous emission factor j3 is not a free 
parameter any more, but comes as a solution to the prob¬ 
lem. We solve the system of coupled differential equa¬ 
tions (1-3) by means of standard ODE solvers imple¬ 
mented in MATLAB. The main advantage of this model 
is that it includes spectral dispersion of scattering, ab¬ 
sorption and gain. 

We choose to simulate typical experimental conditions 
[6]: the sample is a slab 50 /im thick, with 4 = 1.5 ^m, 
and doped with a concentration of ImM of rhodamine 
6G dye whose emission properties are taken from ref. 
[25]. The system is pumped with a 6 ns laser pulse at a 
wavelength of 532 nm. 

Firstly, we consider a conventional RL, such as a poly- 
disperse TiC >2 powder, where i t is not dispersive but 
roughly constant over the gain spectrum. In Fig. 1(a) we 
plot the peak spectral fluence diffusing from the system 
at the air-sample interface together with the full width 
half maximum (FWHM) of the emission spectra as a 
function of the pump energy density. The onset of the 
lasing emission is at P = 0.07 mJ/mm 2 , at higher pump 
intensities the peak fluence increases super-linearly until 
gain saturation is reached; at the same time the emis¬ 
sion width quickly decreases from the broad fluorescence 
emission to a narrow almost constant value at satura¬ 
tion. The threshold of the lasing emission is usually 
identified in analogy to conventional lasing by consid¬ 
ering the change in slope of the pump-peak emission 
intensity relation or by considering the narrowing of the 
lasing emission. Throughout this letter we identify the 
threshold as the half narrowing of the emission spec¬ 
trum, which is a typical experimental parameter more 
accessible than the change of slope [5 ]. The FWHM of 
the emission narrows from the initial ~ 35 nm of the rho¬ 
damine fluorescence spectrum to ~ 0.8 nm of maximum 
lasing narrowing. The fj parameter can be calculated 
from Fig. 1(a) and is /3 ~ 0.04. In Fig. 1(b) the emission 
spectra below and above threshold and for an interme¬ 
diate value are shown. Lasing occurs at the maximum 
of the gain, which is A = 555.5 nm. The width of the 
lasing peak is limited by the losses through the sam¬ 
ple boundaries and by the gain saturation, and, as we 
confirmed, it is not affected by the numerical spectral 
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Fig. 1. (a) Characteristic plot of a random laser: the to¬ 
tal peak intensity (blue squares) and FWHM of the emission 
spectrum (red circles) are plotted for increasing pump energy 
density for a slab system of thickness L = 50 /j,m, transport 
scattering length it = 1.5 pm, and a 1 mM concentration of 
rhodamine 6G providing the gain, excited by a single pulse 
of duration d, = 6 ns at 532 nm. The black dashed line is a 
guide to the eye to highlight the change of regime from flu¬ 
orescence to lasing, (b) The emission spectra (normalised) 
are reported at different pump intensities: the fluorescence 
spectrum (blue) narrows down for increasing pump inten¬ 
sity (red line) down to a 0.8 run FWHM peak (yellow line). 
The corresponding pump intensities are colour coded by the 
hexagons in panel (a). The lasing occurs at the maximum of 
the gain curve (dashed green line). 


discretisation which is set to 0.1 nm. This value dif¬ 
fers from the typical experimental values where the final 
linewidth is typically in the range 5 to 10 nm [6, 26], as 
we are not considering any additional homogeneous and 
inhomogeneous line-broadening effects. 

We consider now the case of a medium characterised 
by resonant Mie scattering. We simulate a system com¬ 
posed of close-packed (n = 1.6, filling fraction / = 0.5) 
dielectric spheres with 1% radial polydispersity. We cal¬ 
culate it by using Mie theory in the approximation of 
independent scatterers [27]. The Mie resonances mod¬ 
ulate 4 as shown by the full lines in Fig. 2(a), where 
scattering is plotted together with the predicted lasing 
peaks for three sphere diameters. Qualitatively, the las¬ 
ing frequency is pulled towards the strongest scattering 
frequency within the gain curve. Compared to what 
is shown in Fig. 1(b), now the lasing frequency can be 
tuned by choosing the Mie resonance. 
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Fig. 2. The effects of resonant scattering on the lasing spec¬ 
trum. (a) The Mie resonances of spherical scatterers (n = 1.6 
in air) of similar size (polydispersty 1%) are exploited to tune 
the lasing emission. The emission wavelength is close to the 
minimum of the transport mean free path while inside the 
gain region (black line), (b) At increasing pump power the 
emission peak wavelength shifts from the fluorescence maxi¬ 
mum to the lasing wavelength stabilizing above threshold. 


Fig. 2(b) shows the narrowing and shifting of the emis¬ 
sion peak while increasing the pump energy. The cross¬ 
ing of the lasing threshold is now evident as a narrowing 
of the emission spectrum while the lasing peak shifts 
away from the maximum of the fluorescence curve. Fre¬ 
quencies close to the maximum of the gain initially dom¬ 
inate, as they receive a larger fraction of the spontaneous 
emission, but those with larger scattering and gain even¬ 
tually prevail above threshold. 

Finally, we study the range of tunability that can be 
achieved by resonant scattering. In Fig. 3(a) we present 
the calculated transport mean free path at A ga i n = 
555.5 nm, which is the wavelength corresponding to the 
maximum of the gain. While for very small scatter¬ 
ers, which are in the Rayleigh regime, the scattering 
increases with the particle size, approaching the Mie 
regime this trend is reversed. The minimum scatter¬ 
ing length is achieved at r < A/2, which corresponds to 
the onset of the first Mie resonances, and subsequently 
increases with a resonant behaviour, while the res¬ 
onances get closer and closer. As previously shown in 
Fig. 2, this modulation in the scattering induces a tun¬ 
ing of the emission wavelength which is now explicit in 
Fig. 3(b) where we plot the lasing frequency versus the 
particle reduced radius r/X ga i n . While in the Rayleigh 
regime the lasing emission follows the maximum of the 
gain, for r > A/2 the emission starts to red-shift follow¬ 


Fig. 3. Tuning the lasing emission wavelength via reso¬ 
nant scattering, (a) Transport scattering mean free path 
It computed by Mie scattering (spheres of refractive index 
n = 1.6 in air, filling ratio / = 0.5 and 1% polydispersity) 
at Xmax = 555.5 nm. (b) Lasing emission wavelength for 
resonant scattering. When the particle size is increased from 
non-resonant Rayleigh scattering (r << A) to Mie scattering, 
the lasing emission wavelength is tuned to follow the scatter¬ 
ing resonance in a range of 25 nm inside the gain curve (gray 
area). The red lines are further plotted in Fig. 4(a). 


ing the Mie resonances. The wavelength jumps corre¬ 
spond to the appearance of a more favorable resonance 
which pulls the lasing wavelength in the blue part of the 
gain region. The tunability range for the simulated sys¬ 
tem is highlighted by the area in grey and it is roughly 
25 nm. 

Fig. 4(a) is a zoom in of Fig. 3(b) and shows in detail 
the lasing frequency evolution around r = A. The curve 
consists of three different continuos branches separated 
by two jumps. The numerical results of the model are 
compared to analytical diffusion calculations which pre¬ 
dict lasing occurring at the wavelength which minimises 
the critical length L cr [6, 7], the length above which the 
gain exceeds the losses. In a slab geometry, and under 
the approximation of uniform gain, the critical length is 
L cr = TTy/£ t l g / 3, where l g = The minimum 

of L cr and the numerical results are in good agreement, 
with significant differences only close to the wavelength 
jumps where the simplified analytical model fails to ac¬ 
count for the competition for the gain. We attribute 
this to the increased mode competition between two dif¬ 
ferent Mie resonances. Mode competition is strong for 
particle sizes around r = 571 nm and 595.5 nm, and 
this is reflected by an increased threshold. For instance, 
the threshold does not always decreases monotonically 
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Fig. 4. Details of the emission wavelength and lasing thresh¬ 
old. (a) Pulled by the Mie resonances, the lasing emission 
redshifts until the next Mie resonance enters the gain spectral 
region. This competition is shown by the tooth-saw profile. 
The numerical values (squares) are compared to the analyti¬ 
cal prediction (orange line) obtained by the minimum of the 
critical length L cr = Tv^/itlg/3. (b) The threshold, defined 
as the point of half emission narrowing, is quite constant in 
the tuning range but close to the frequency jumps, when the 
competition between the two lasing peaks strongly increases 
it. 

when the lasing emission approaches A ga in- For the 
other lasing frequencies the threshold is quite constant 
in the whole region of tunability with values around 
0.1 mJ/mm 2 . In addition, we confirm (not shown here) 
that for particles size r = 571.9 nm and r = 595.4 nm 
the spectrum develops two competing peaks of which 
one finally prevails for large enough pump energy. 

In conclusion, a dispersive diffusive gain model is ca¬ 
pable to simulate RL action when driven by gain or scat¬ 
tering resonances. For resonant scattering we calculated 
a 25 nm emission tuning range, with minimal threshold 
increase. Our dispersive model does not require any as¬ 
sumption about the /3 parameter and can predict mode 
competition and a non-trivial threshold dependence. As 
a further extension, the model can include arbitrary ab¬ 
sorption curves and dispersion also for the energy veloc¬ 
ity, and can be extended to the full three-dimensional 
case at the expense of increased computing time. 
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